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ABSTRACT 

High resolution 2D hydrodynamical simulations describing the evolution of the hot ISM in 
axisymmetric two-component models of early-type galaxies well reproduced the observed 
trends of the X-ray luminosity (Lx) and temperature (Tx) with galaxy shape and rotation, 
however they also revealed the formation of an exceedingly massive cooled gas disc in rotating 
systems. In a follow-up of this study, here we investigate the effects of star formation in the 
disc, including the consequent injection of mass, momentum and energy in the pre-existing 
interstellar medium. It is found that subsequent generations of stars originate one after the 
other in the equatorial region; the mean age of the new stars is > 5 Gyr, and the adopted recipe 
for star formation can reproduce the empirical Kennicutt-Schmidt relation. The results of the 
previous investigation without star formation, concerning Lx and Tx of the hot gas, and their 
trends with galactic shape and rotation, are confirmed. At the same time, the consumption of 
most of the cold gas disc into new stars leads to more realistic final systems, whose cold gas 
mass and star formation rate agree well with those observed in the local universe. In particular, 
our models could explain the observation of kinematically aligned gas in massive, fast-rotating 
early-type galaxies. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: ISM - galaxies: kinematics and 
dynamics - ISM: evolution - X-rays: galaxies - X-rays: ISM 


1 INTRODUCTION 

Early-Type galaxies (ETGs) are embedded in hot gaseous haloes 
produced mainly by stellar winds, and heated to X-ray temperatures 
by Type la supernovae (SNIa) explosions and by the thermalization 
of stellar motions jFahbiano| 1 989| |Q’Sullivan et al.|2001[|Ciotti| 
|et al.|199i~| |David et al.|1991} . The thermalization is due to the 
interaction between the stellar and SNIa ejecta and the pre-existing 
hot ISM (e.g., |Mathews|l989j[Parriott & Bregman|2008[ l. Recent 
high resolution 2D hydrodynamical simulations of hot gas flows 
j Negri et al.|2014a||Negri et al.|2014b| hereafter N14) showed that 
the presence of ordered rotation in the stellar component can alter 
significantly the ISM evolution with respect to that shown by fully 
velocity dispersion supported systems of same total mass and mass 
distribution. Firstly, it is found that the rotation field of the ISM in 
rotating galaxies is very similar to that of the stars, with a consequent 
negligible heating contribution from thermalization of the ordered 
motions. Secondly, conservation of angular momentum in the ISM of 
rotating galaxies results in the formation of a centrifugally supported 
cold equatorial disc, with the consequent reduction of both the X-ray 
luminosity Lx and temperature Tx of the hot ISM. These results 
compared well with observations, which show a dependence of L x 
and T x on the galactic shape and internal dynamics: L x is observed 
to be high only in round and slowly rotating galaxies, and is limited 
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to lower values for flatter, fast rotating ones (|Eskridge et al.|1995| 


Pellegrini et al.| 1 997]|Sarzi et al.|2010||Pellegrini|2012[ Li et al.| 


2011 ; |Sarzi et al.|2013| ); Tx of slowly rotating systems is consistent 

just with the thermalization of the stellar random kinetic energy, 
estimated from cr e (the stellar velocity dispersion averaged within 
one effective radius R e ), while fast rotating systems show Tx values 
below 0.4 keV, and not scaling with cr c jSarzi et al.|2013| see also 
|Pellegrini|201 l|[Posacki et al.|2013| >. 

A major outcome of the previous simulations is that cold mate¬ 
rial can be accumulated in considerable amounts, during the lifetime 
of rotating galaxies. The cold gas typically settles in the equatorial 
plane, where it forms an extended disc (of 0.5-10 kpc radius), that 
can be as massive as ^ 10 10 Mq, in the largest galaxies. This result 
brings in a few important questions: are these cold discs observed? 
do they become, as seems natural, a site for star formation (here¬ 
after SF)? is SF observed? and what is the impact of SF in the disc 
on Lx and L x ? In recent years evidence has been accumulating 
that ETGs host significant quantities of cold gas, in the form of 
atomic and molecular hydrogen (Morganti et al. 2006, Combes et al. 


di Serego Alighieri et al.|2007| Grossi et al.|2009| Young 


2007 


e t al.|2011) ; approximately 50% of massive ETGs (of stellar mass 

M, > 10 rtl Mg) contain 10 7 to 10 9 Mg of HI and/or H 2 {Young et al. 

[20T4l [Serra et al.|2012|[2014) |Davis et al.|201 l||20l3] T A large, 

systematic investigation of the ATLAS 30 sample of ETGs with the 
Westerbork Synthesis Radio Telescope found that ^ 40% of galaxies 
outside Virgo, and ^ 10% of galaxies inside it, are detected in HI, 
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with M H i > 10 7 Mq; the majority (2/3) of the detections consists 
of settled configurations, where the cold gas is in discs or rings. 
Small discs (size of a few kpc), confined within the stellar body, 
share the same kinematics of the stars; large discs (up to 5 x 10 9 
Mq) extend to tens of kpc, and in half of the cases are kinemat¬ 
ically misaligned with the stars ( |Serra et al.|20 1 2\ . In particular, 
fast-rotating Virgo galaxies have kinematically aligned gas, and the 
most massive (M, > 8 x 10 10 Mq) fast-rotating ETGs always have 
kinematically aligned gas, independent of environment; this align¬ 
ment leads to hypothesize that the gas has been internally generated 
{Davis et al.|20 1 1 [|20 13| . The HI discs/rings around slowly-rotating 
ETGs, instead, are usually not fully settled, which suggests an exter¬ 
nal origin (in mergers, or in accretion from satellite galaxies, or from 
the intergalactic medium); an external origin was considered likely 
also for those ETGs showing stellar/gas kinematic misalignment. 
Interestingly, while HI is more ubiquitous, molecular gas is detected 


only in fast rotators across the entire ATLAS 30 sample Young et al. 
[20TT|[20T4l|Davis et al.|20TTl[20l3l >. 


Besides giving indication about its origin through its morphol¬ 
ogy and kinematics, the observed cold gas seems also to provide 
material for SF. Low level SF activity is present in ~ one-third of 
ETGs {Yi et al.|2005l|Suh et al.|20l0l|Ko et al.|2014| >; |Sarzi et al.| 
( |2006[ > showed ongoing SF signatures in the optical spectra of at 
least ~ 10% of nearby ETGs. In the ATLAS 30 sample, galaxies with 
HI within ~ 1 R e exhibit ongoing SF in ~ 70% of the cases, ~ 5 
times more frequently than galaxies without HI {Serra et al.|20 12{ . 
Interestingly, as for molecular gas, some degree of SF and young 
stellar populations are detected only in fast rotators, in the ATLAS 30 
sample {Kuntschner et al.[2010[|Sarzi et al.|2013| >. Integral-held 
spectroscopy showed that ETGs host frequently a rotating stellar 
component younger and more metal rich than the bulge ( |Krajnovic| 
|et al.|2008} . The presence of this component, and the occurrence 
of SF, imply an important role for the cold gas during the evolution 
of ETGs {Khochfar et al. |20 1 l||Naab et al.|20l4[|Cappellari et al.| 
|20l3l >. 


Finally, the cold gas content of ETGs is also an important pre¬ 
diction of ACDM hydrodynamic al simulations of galaxy formation, 
that include cold gas evolution (e.g., |Martig et al.|2013{ , and accre¬ 
tion via various processes during the secular evolution of galaxies 
{Oser et al.|2010| see also |Lagos et al,|2014||Dubois et al,|20 13 j i . 

In conclusion, the cold gas has become a tool to gain insight 
into recent (and less recent) galaxy evolution. In order to correctly in¬ 
terpret the variety of observational results, and to use them properly 
as constraints for different scenarios for the origin of the structure 
and SF history of ETGs, it is crucial to establish what is the relative 
importance of the various gas production processes (internal and ex¬ 
ternal to galaxies), gas depletion ones (AGN and SF-driven outflows, 
environmental stripping), and gas consumption ones (SF). Cold gas 
is usually thought to come from accretion from the surrounding 
medium or satellites, as well as from gas-rich mergers; in the cases 
of the giant central-dominant galaxies in groups or clusters it can 
also come from cooling of hot gas (Edge et al. 2010, McDonald et al. 
2011). The numerical investigation of N14 showed an additional 
internal contribution of cold gas, coming from the evolution of the 
passive stellar population, that can be substantial, or even too large 
with respect to observed values. A natural sequel of the N14 work 
should address the questions of whether this gas can possibly lead 
to SF, when SF takes place, and what is the fate of the cold discs, 
whether they are consumed or they are continuously replenished by 
cooling hot gas. In this paper we add SF to the simulations of N14. 
and we explore the ISM evolution including the removal of cold gas, 
and the injection of mass, momentum and energy appropriate for the 


newly (and continuously) forming stellar population. In this way we 
aim at establishing whether i) the N14 results for the general trends 
of the hot gas properties with galaxy shape and stellar kinematics 
still hold; ii) the formation of stars can reduce the amount of cold 
gas in the simulations, thus bringing it more in agreement with ob¬ 
served values; iii) a significant channel for SF, previously neglected, 
should be taken into consideration for rotating systems, and whether 
this can account for the low level SF activity currently seen to be 
ongoing. 

This paper is organized as follows. In Section 2 we describe the 
main ingredients of the simulations, such as the galaxy models and 
the input physics. In Section 3 we present and discuss the results of 
the simulations. In Section 4 we summarize the conclusions. 


2 THE SIMULATIONS 

N14 performed a large set of 2D hydrodynamical simulations with 
the ZEUS MP2 code to fully explore the large parameter space of 
realistic (axisymmetric) galaxy models, characterized by different 
stellar mass, intrinsic flattening, distribution of dark matter, and 
internal kinematics. The galaxy flattening was either fully supported 
by ordered rotation, originating the set of models that are isotropic 
rotators, or by tangential anisotropy, originating the set of fully 
velocity dispersion supported models. These two extreme configu¬ 
rations were built adopting the Satoh decomposition, respectively 
with Satoh parameter k = I and k = 0. The galaxy models were 
tailored to reproduce the observed properties and scaling laws of 
early type galaxies (see also [Posacki et al.p0l3] l. In this work we 
perform hydrodynamical simulations for a representative subset of 
rotating models already investigated by N14 (Sect. 2.1), including 
SF in the simulations (as described in Sect. 2.2 below), but keep¬ 
ing the code (numerical set-up, grid properties) in all equal to that 
used by N14 (see N14 for details). Also, a less extreme value of 
k = 0.1 is explored. A logarithmically spaced numerical mesh ( R , z) 
of 960x480 gridpoints is employed, with a resolution of 90 pc in 
the first 10 kpc from the centre. 


2.1 The galaxy models 

N14 built axisymmetric two-component galaxy models where the 
stellar component has two different intrinsic flattening, correspond¬ 
ing to the E4 and E7 shapes, while the dark matter halo is spherical. 
The luminous matter is described by the deprojection {Mellier &| 
Mathez 1987) of the[de Vaucouleurs (1948) law, generalized for 
ellipsoidal axisymmetric distributions; the dark matter profile is the 
( Navarro et al. |1997[ l one, with the dark mass M h amounting at =: 20 
times the total stellar mass M,. For any fixed galaxy mass and shape 
(E4 or E7), N14 built two models: the first one, called “FO-built”, 
when seen face-on has the same R c of the spherical E0 counterpart, 
thus its stellar mass distribution becomes more and more concen¬ 
trated than in the E0 model, as it gets flatter; the second one, called 
“EO-built”, when seen edge-on has the same circularized R s of the 
E0 counterpart, which makes its stellar mass distribution to expand 
with increasing flattening. 

In this work we re-simulate a few flat (E4 and E7) rotating 
models (isotropic rotators, with k = 1) of N14 including SF in 
them. In order to explore the effects of SF at the high and low 
ends of the galaxy mass range explored by N14, we choose four 
flat models with luminosity-weighted stellar velocity dispersion 
within /? e/8 of cr e8 = 300 km s~* for the parent spherical model 
(and we call these high-mass models, “HM” models), and four with 
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Table 1 . Main properties of the galaxy models. 


Name 

Lb 

Re 

M, 

M h 

0"e8 

/dm 

c 


(IO^Lq) 

(kpc) 

(10“M Q ) 

(10 U M G ) 

(km s _I ) 



a) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

EO4 2 s 00 

0.26 

4.09 

1.25 

25 

166 

0.63 

37 

E07[ s °° 

0.26 

4.09 

1.25 

25 

124 

0.66 

37 

F04 2 °° 

0.26 

4.09 

1.25 

25 

178 

0.59 

37 

F07[ s °° 

0.26 

4.09 

1.25 

25 

150 

0.57 

37 

EO4jJ 0 

0.62 

7.04 

3.35 

67 

207 

0.62 

28 

EO7jJ 0 

0.62 

7.04 

3.35 

67 

154 

0.66 

28 

EO4[ s 00 

1.32 

11.79 

7.80 

160 

248 

0.64 

22 

E07[ s °° 

1.32 

11.79 

7.80 

160 

185 

0.68 

22 

FO4 3 ° 0 

1.32 

11.79 

7.80 

160 

266 

0.60 

22 

F07 3 °° 

1.32 

11.79 

7.80 

160 

224 

0.59 

22 


Notes. (1) Model name: “EO” or ‘TO” indicate the procedure (EO-building or FO-building) applied to the spherical (EO) model to obtain the shape indicated by 
the number (E4 or E7 shapes); the superscript indicates cr e 8 of the corresponding EO model; the subscript “IS” indicates the kinematical configuration of the 
isotropic rotator (for example, F04“® 0 is an E4 isotropic rotator, FO-built from a spherical model with cr c 8 = 200 km s 1 ). (2) Luminosities in the B band. (3) 
Effective radius (for a FO view for FO-built models, and an EO view for EO-built models). (4) Total stellar mass of the original (old) stellar population. The first 
four models are LM models, the next two are intermediate mass models, the last four are HM models. (5) Total dark matter mass. (6) Luminosity-weighted 
average of the stellar velocity dispersion within a circular aperture of radius R e /$] f° r non-spherical models, <x e 8 is the edge-on viewed value. (7) Dark matter 
fraction enclosed within a sphere of radius R e . (8) Concentration parameter of the dark matter profile. 

Note: the properties listed above, for the E04 250 and E07 250 models with k = 0.1. are not reported, since they are equal to those of the E04[J° and E07[J° ones, 
except for cr e 8. which is respectively 223 km s -1 for the E04^® l model, and 184 km s _1 for the E07|^ t model. 


<r e8 = 200 km s” 1 for the E0 counterpart (low-mass models, “LM” 
models). In N14, the first set was found to host inflows, with the 
creation of a massive, centrifugally supported cold gaseous disc, 
while the second set was found in a global wind, or close to the 
transition to it. In a global wind the gas has very low density and 
positive velocity (directed outwards) through most of the galaxy 
(e.g., Mathews & Baker 1971 [[ In addition, we re-simulated two 
intermediate-mass models (again E4 and E7), with <x e8 = 250 km 
s” 1 for the E0 counterpart. For these two models, we also built 
a moderately rotating stellar kinematical configuration (k = 0.1), 
not explored by N14. The main structural properties of the ten re¬ 
simulated models, identical to those presented in N14, are listed in 
Tab. 1; the two new models with k = 0.1 differ from the intermediate- 
mass models only in the tr c8 value, which is given in the notes to 
the table. 


2.2 Star formation in the code 

We employed two different schemes for SF, a passive (pSF) one 
and an active (aSF) one. In the pSF, only cold gas removal from the 
numerical mesh is allowed; in the aSF, we also consider the injection 
of mass, momentum and energy from the newly forming stellar 
population, for simplicity limiting in this work to the evolution of 
stars more massive than 8 Mq (ending with SNII explosions); note 
however that for reasonable stellar initial mass functions (IMF) these 
massive stars are major contributors to the total mass return rate. We 


1 In fact the phase of the gas flow (that can range from a wind to an inflow) 
is basically determined by the different relative importance of SN heating 
in galaxies of different mass, as already thoroughly discussed injCiotti et al.| 
[1991] 


describe below the main inputs and sinks of mass and energy for the 
gas flow, due to both the original stellar population (of total mass 
M.) and the newly forming stars, and the corresponding equations 
of hydrodynamics solved by the code. 


2.2.1 The mass injection and sink terms 

The mass inputs are stellar winds and SNIa’s ejecta produced during 
the passive evolution of the original stellar population of the galaxy 
(at a rate per unit volume respectively of p, and p Ia , for which we 
adopt the standard recipes coming from the stellar evolution theory; 
e.g.. N14), and the Type II supernovae ejecta produced by the newly 
born stellar population (at a rate of p n , calculated as detailed below). 

The mass sink is due to SF subtracting gas from the grid, at an 
adopted rate per unit volume of: 

77sfP 

Psf = --, f S F = max(f coo i, f d yn), (1) 

‘SF 

where p is the gas density, and t] SF is the SF efficiency, for which we 
adopt two values of t/ sf = 0.01 and 0.1. The rate p SF depends on the 
maximum between the cooling timescale f coo i = E/X (where E is the 
ISM internal energy density, and X is the ISM bolometric luminosity 
per unit volume), and the dynamical timescale f dyn = f'inl'ilGp 
(see also |Ciotti & Qstriker|2007| . 

When SF takes place in a computational cell, the ISM is re¬ 
moved and an equal mass in stars AM, appears at the same place. 
These new stars are assumed to form with a Salpeter IMF, thus for a 
given AM,, the number of stars having a mass greater than 8 Mq, 
that will explode as SNII, is Nu - 7 x 10” 3 AM,(Mq). These new 
stars in turn inject mass into the ISM. By integration of the mass dif¬ 
ference between the mass of the progenitor and that of the remnant, 
for the Salpeter IMF, one finds that the mass injected by SNITs is 
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=: 0.2AM,. The final SNII mass source term p n comes from consid¬ 
ering that a given SF episode generates SNII’s that inject mass at 
a rate exponentially declining on a timescale tu = 2 x 10 7 yr, and 
that at a certain time during the evolution of a SF episode, another 
episode may take place, forming younger SNIFs that in turn eject 
mass into the ISM. The whole process can be formalized as: 

dpii _ _Pu + 0-2psF 
d? Tu Tu 

( |Ciotti & Qstriker|2007j i. The equation above is very useful in nu¬ 
merical works, as it allows to compute the mass return of the new 
stars formed at each timestep without storing the whole SF history at 
each gridpoint, but only the current value and the value of the mass 
return at the previous timestep. For more complex recipes, built on 
the scheme above, see|Calura et al.||2014}. 


Finally, an energy and momentum sink associated with SF are 
present; they are respectively written as: 


Esf - 


77sf£ 


m sf 


risFin 


( 8 ) 


tsF t S f 

where E and m are the internal energy and momentum density of 
the ISM. 


2.2.3 The hydrodynamical equations 

The hydrodynamical equations are the same as in N14, with the 
addition of all the source and sink terms related with the SF process 
described above: 

3p 

—— l- V • ( pu ) = pi a + p» + pn - Psf, 
ot 


2.2.2 The energy injection and sink terms 


Energy is injected into the ISM by the thermalization of the kinetic 
energy of SNIa and SNII explosions, and by the thermalization of 
random and streaming stellar motions for the stellar winds. The rate 
of SNIa’s explosions is the same entering p ia defined in Sect. |2.2.T| 
The energy input rate per unit volume from the original (old) stellar 
population is then: 


Eon - E\.. + 


Pia + P* 


I ||2 o 

v v e„ - h + Tr(o^) 


(3) 


where £j a = Pi a i?sNia£sN/l-4 Mq, $ SNIa is the thermalization effi- 
ciency (for which we ado pt the value of 0.85, see |Thomton et al.| 


1998 


Tang & Wang 2005 1 , E s N = 10 51 erg, and 1.4 Mq is the mass 


that is ejected by one SNIa event; v v is the stellar streaming velocity 
field, u is the velocity of the ambient gas, tr 2 is the stellar velocity 
dispersion tensor of the stars. 

The energy input rate £ new from the newly bom stellar popula¬ 
tion is derived as follows. A mass of newly formed stars AM, injects 
energy through SNII with an efficiency 


£ll 


AjrEsN 
A M,c 2 


=* 3.9 x 1(T 6 . 


(4) 


Consistently with eq. [ 2 ] the SNII energy injection rate En due to the 
thermalization of the ejecta is given by: 


dEn _ £j] £iiC - Psf#snii 

1. ~ ' ’ W/ 

at r„ Tn 

where we take the fiducial value of i?snii = #SNia/5 = 0.17 to 
account for the lower thermalization efficiency of SNII exploding in 
a cold and dense medium. We assume that the new stars inherit the 
kinematical configuration of the original stellar component in the 
place where they are born (i.e. the velocity dispersion and rotation of 
the new stars are the same as those of the original stellar distribution 
in the same place). This allows us to treat the heating terms described 
in eq. (7) below by using the same properties of the old stars. Note 
however that this choice is not unreasonable, because the bulk of 
star formation takes place in the cold and rotationally supported disc, 
whose rotational properties are very similar to those of the stars. 
Thus, the total heating rate per unit volume due to Type II SNe is 


Enew = -Ell + Y - m||“ + Tr(£T 2 )J . (6) 

The total energy injection due to the old and new stellar populations 
is then E = £ 0 i t | + E ntvl , that is: 


E = 


-Fla + £jl + 


Pia + P, + PlI 
2 



(V) 



+ p(u ■ V)« = -Vp - pVO tot + 


+ (Pia +P. +Pn)OA - »)> 


- t- V ■ (Eu) = -pV • it - _£ + E - E s p, 

ot 

where p, u, E, p, O tot , and £. are respectively the ISM mass den¬ 
sity, velocity, internal energy density, pressure, total gravitational 
potential, and bolometric luminosity per unit volume. The gas is 
assumed to be an ideal monoatomic fully ionized plasma, so that 
p = (y - 1 )E, where y = 5/3 is the adiabatic index. The chemical 
composition is fixed to solar (p = 0.62), and the gas self-gravity is 
neglected. 

In the pSF scheme, the source terms pn and En are zero, thus 
they either do not enter the hydrodynamical equations above and 
they do not contribute to E. In the more realistic aSF, p n and En are 
non-zero. Most galaxy models are simulated with pSF and aSF; for 
each of the two schemes, we adopted two values of the SF efficiency 
in eqs.[l]and[8] 77 SF = 10 _1 and p S F = 10~ 2 (see Sect. 3 below). 

We tested that the code conserves the total mass and energy; 
the total mass at any time t is M gas + M esc = M m j + MP. - M" ew , 
where M gas is the gas mass within the simulation box, M esc is the 
cumulative escaped mass out of the simulation box, until that time; 
M inj , Mf and Mf w are the cumulative masses respectively injected 
by the passively evolving stellar population, injected by SNII events 
due to the newly born stellar population, and of the new stars. The 
values of these masses at the end of the simulations are given in 
Tabs. |2|5] Finally, the radiative cooling is implemented by adopt¬ 
ing a modified version of the cooling law of jSazonov et al.| ( j2005| , 
with a lower limit for the ISM temperature of T = 10 4 K. The 
total X-ray emission in the 0.3-8 keV band (Lx), and the emission 
weighted temperature in the same band (7x). are calculated as vol¬ 
ume integrals over the whole computational grid, using as weight 
the emissivity in the 0.3-8 keV band of a hot, collisionally ionized 
plasma (see N14 for more details). 


3 RESULTS 

Each simulation starts with the galaxy empty of gas, when the age 
of the original stellar population is 2 Gyr, and the ISM evolution 
is followed for the next 11 Gyr. Tables [2]5] list for all models the 
main quantities of interest at the end of the simulations; the tables 
also list the same quantities for the corresponding models without 
SF taken from N14 (hereafter “noSF” models). First (Sects. HU we 
present the main results focussing on the two sets of re-simulated 


© 2015 RAS, MNRAS POQ.fTp)] 





















X-ray haloes and star formation 5 



150 200 250 300 150 200 250 300 150 200 250 300 

cr e 8 (km s -1 ) ^8 (kms' 1 ) cr e8 (kms -1 ) 

Figure 1. Top panels: ISM X-ray luminosity Lx in the 0.3-8 keV band at 13 Gyr as a function of cr e 8, for the selection of ten models from N14 re-simulated 
here (Tab.^, plus two models with cr e g = 250 km s -1 for the E0 counterpart, and k = 0.1 (shown with thicker ellipses). The three panels (from left to right) refer 
to the same models without SF (from N14), with passive SF, and with active SF (with 77 SF = 0.01). Bottom panels: the same for the X-ray emission weighted 
temperature Tx in the 0.3-8 keV band at 13 Gyr. All Lx and Tx values are given in Tabs. |2|5] See Sect. |3.1| for more details. 



0.2 0.4 0.6 0.8 

T x (keV) 


Figure 2. Lx as a function of Tx for aSF models listed in Tables [I|5l (open 
ellipses, as in Fig. [TV compared with Chandra data for the hot gas of ETGs 
taken from Boroson et al. < 2011 1 and Kim & Fabbiano <20151 (grey squares 
with error bars). From the latter two samples we selected only ETGs of 
morphological type later than E3, for a proper comparison with our models. 


models with largely differing values for the galaxy mass (i.e., the 
two sets with E0 counterpart of <x e8 = 200 and 300 km s -1 , the LM 
and HM models). In Sect. |3.2| we concentrate on the intermediate 
mass models (run only for the aSF), where this time k is equal to 1 or 
lower (k = 0.1). In Sects. [33| and |3.4| we discuss the consumption of 
cold gas mass and the formation of new stars. Finally, we investigate 
the relationship between the adopted recipe for SF (eq.[7} and the 
Kennicutt-Schmidt relation (Sect. [33}. 

3.1 Comparison between noSF and SF models 

In a short summary, for noSF models, N14 found that rotation in 
LM galaxies favours the establishment of global winds, with the 
consequent reduction of Lx ; in medium-to-high mass galaxies the 
conservation of angular momentum lowers the hot gas density in the 
central galactic region, leading again to a reduction of L x , and also 
of Tx (because the external and colder regions weight more in the 
computation of T x ). In LM galaxies, instead, T x can become higher 
if rotation triggers a wind, due to the decrease of the ISM density, 
and the additional heating due to the high meridional velocities of 
the escaping material (see eq.[7}. 

The first important result here is that when SF is added (both 
in the pSF and aSF modalities), the hot ISM evolution remains 
substantially similar to that found by N14 without SF. for models 
of the same mass. As a consequence of the insensitivity of the 
general behaviour to the addition of SF, the values of L x and T x 
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Figure 3. Meridional sections of the heating over cooling time ratio for the LM E07^° model, at a selection of representative times (indicated as galaxy ages 
in the top right of each panel). The heating time is ?heat = E/E (where E is the source term given in eq.[7}; the cooling time is defined in Sect. 2.2.1. From 
left to right, each column refers to the noSF, pSF and aSF (with t/sf = 0.1) models, respectively. Arrows describe the meridional velocity field; their lenght is 
proportional to the modulus of the gas velocity in the (z, R ) plane, according to a scale shown in the bottom left comer. See Sect. [XT] lor more details. 
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Figure 4. Meridional sections of the heating over cooling time ratio for the HM E07^° model. From top to bottom, each row refers to the noSF, pSF and aSF 
(with r/sF = 0.1) models, respectively. Arrows describe the velocity field, as in Fig. [3] Note the large cold disc, that does not completely disappear with SF. See 
Sect. l3Tl for more details. 


for the models with SF. and their trends with the main galactic 
structural properties, are on average very similar to those found by 
N14. Therefore, the conclusions of N14 about the importance of 
shape and rotation in determining L x and T x , also as a function 
of galaxy mass, are confirmed. Figure [T] demonstrates this result 
by plotting the final L x and T x values for the same rotating galaxy 
models evolved with noSF and with SF (with rj S p = 0.01): the 
distribution of the points in the noSF and SF panels is very similar 
(and the analogous figure with the = 0.1 models produces 
the same conclusions). In particular, L x of HM models is only 
marginally sensitive to the presence of SF, while in LM models 
SF can introduce variations of L x ; however, these keep within the 
(large) range of L x values already found in N14 without SF. Thus, 
SF adds another cause of spread for L x , at lower galaxy masses. 
Figurc[2]shows the resulting L x vs. T x for the aSF models, compared 


with the same quantities derived recently from Chandra data for 
the hot gas of ETGs ( |Boroson et al.|2011| and |Kim & Fabbiano| 
|2015[ >. From the latter two samples we have selected only ETGs of 
morphological type later than E3, for a proper comparison with our 
models. Overall, the simulation results agree well with the observed 
X-ray properties of real ETGs (except perhaps for the most X-ray 
luminous models that may be slightly hotter than observed). 

The reason for variations in L x linked with SF in LM models 
is given by anticipations or delays in the most conspicuous features 
of their typical flow behaviour: N14 found (and the present simula¬ 
tions confirm) that the ISM in rotating models experiences periodic 
cooling episodes, where the gas injected by the stellar population 
accumulates until the radiative losses become catastrophic. In these 
episodes, the ISM quickly cools, emitting a large part of its internal 
energy as radiation; peaks in L x and throats in T x are then produced, 
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Figure 5. Time evolution of the X-ray luminosity Lx and X-ray emission weighted temperature 7x for the E07gr (top panels), F07jg (middle panels), 
(bottom panels) models. The black, green and red lines refer to the noSF model, and to the passi ve an d active star formation schemes, respectively, with 

for more details. 


ECW ;™ 0 


>7SF 


: 0.1. The F07j®° model has an evolution almost identical to that of E07j®° shown here. See Sect. 


3.1 


as apparent in Fig. 5. The tenuous hot atmosphere left after a major 
cooling is then replenished by the new mass injection from the stellar 
population. In correspondence of the formation of these significant 
amounts of cold material (with short cooling times) star formation 
is enhanced, and peaks in pace with L x , as discussed in Sect. 3.4. 
On the contrary, very little variations of the flow due to SF are seen 


in HM galaxies. Figures [3] [4] and [5] give particular examples of 
the general similarity in the hydrodynamical evolution of noSF and 
SF models, by plotting in parallel the ISM evolution for the noSF. 
pSF and aSF cases; at the same time, the figures point out how this 
similarity is very close for HM models (see Fig. [4] for the HM model 
E07j^°), while some differences can be present in LM ones (see 
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Fig. 0 for the LM model E07™). Notably, the hot ISM evolution 
in the HM models is practically identical with and without SF (see 
also Fig. [5}: a massive cold disc forms early, and afterwards, even 
when it is mostly consumed by SF (see Sect. |3.3| below), around the 
disc the hot ISM evolution proceeds almost unaltered by the energy 
and mass injection from SNII’s. Less massive models, instead, are 
sensitive to the inclusion of SF, similarly to what found in previous 
studies where they showed to be very sensitive to any change in the 
galaxy properties (e.g., the mass distribution, the stellar population 
inputs, the stellar kinematics; Ciotti et al .| 1991 1 . N14). As shown 
by Figs.[3]and[5] for the E07 2 ™* galaxy, at around 4 Gyr the noSF 
and the pSF models have already past their major cooling episode 
(see the peak in L x in Fig. [5J followed by the onset of a wind, and 
the equatorial outflow becomes stronger and stronger thereafter (see 
how the purple regions become more and more extended in Fig. [3] 
starting from 4.6 Gyr, and the sharp and steady decline in L x in 
Fig. The aSF model, instead, has its major cooling catastrophe 
much later, as shown by the prominent dark green region in Fig. [3] 
at = 8 Gyr, which emits a conspicuous amount of radiation in the 
X-ray band, thus prompting the corresponding peak in L x in Fig. [5] 
This delay in the major cooling episode, preceding the onset of a 
wind, is due to the newly formed stars that inject mass, and most 
importantly heat, in the ISM; thanks to this heat the galaxy can 
remain hot gas rich until = 8 Gyr, when finally the major cooling 
takes place. pSF, where mass is not injected instead, corresponds to 
an evolution very similar to that without SF. 

A complementary illustration of the different behaviour be¬ 
tween the HM and the LM galaxies is given by Fig. [5] that shows 
the Lx and T x evolution for the same models of Figs. [3] and [4] 
plus another set of three models representative of the LM class, the 
F07™ ones. In LM models, a number of cooling episodes (peaks 
in Lx) typically take place, possibly terminated by a major one that 
is followed by the onset of a wind and the clearing of the ISM from 
the galaxy (as in N14). The galaxy mass distribution in the FO7jj 0 
models (middle panels of Fig. [5j. is more concentrated than in the 
E07jg° ones, thus the last major peak in L x is delayed with respect 
to what happens for the E07™ models. Thus, in the F07™ galaxy, 
aSF produces a delay in the major cooling episode that is even longer 
than for the E07™ aSF model, to the point that it does not take 
place within the present epoch. pSF instead produces its anticipation 
with respect to what shown in the noSF case (Fig#, due to pSF 
subtracting gas and leaving a lower density region, radiating less 
and more easy to push out of the galaxy. 

A final, general result is also that the evolution of a model with 
SF is more and more different from that of its corresponding noSF 
model, when increasing r/ SF . 

For illustration purposes. Fig. [6] shows a meridional section 
of the density distribution of the newly formed stars at the end of 
the simulations, for the models in Fig. [5] One can notice the more 
extended and massive disc in HM models, where it is also quite 
independent of the pSF or aSF scheme. On the contrary, in LM 
models, the disc is much less extended, and it can become larger 
(for the E07™ model) or smaller (for the FO7j| 0 model) when 
switching from the pSF to the aSF scheme; this is another evidence 
of the high sensitivity of the flow evolution to any change in the 
input parameters, at low galaxy masses. 

3.2 Intermediate mass ETGs 

Here we present the results for two representative models, of E4 
and E7 shapes, with intermediate galaxy mass (cr e g = 250 km s -1 
for the EO counterpart, see Tab. |TJ, that have been run with the 


same kinematic configuration of the LM and HM models presented 
previously (k = I, isotropic rotator), and also with a lower rotation 
(k = 0.1). For these, only aSF has been considered, with ;/sf = 0.01 
and 77 s F = 0.1. The final values for the most relevant quantities of 
these eight models are presented in Tabs.[4]an d0 

The E4 models with k = 1 show the same smooth evolution 
of the HM models, i.e., no or very small oscillations in L x and T x 
(see the dashed lines in Fig. [7] similar to those in Fig. [5j; in fact, 
the evolution of the hot gas is similarly not affected by SF. The 
E7 models with k = 1, instead, show oscillations in L x and T x , 
similar to those of the £07™ and F07™ models, but fewer and 
less pronounced (see the solid lines in Fig. [7] to be compared with 
those in Fig. [5}. In the aSF model the final L x is larger than for the 
noSF case, as found for the F07™ models (Sect. 3.1). 

In models with rotational velocities reduced by a factor of ten 
(k = 0.1), L x and T x increase largely, as expected from the results 
of N14 showing a larger L x and T x in non-rotating galaxies (k = 0). 
Again the E4 models show no oscillations in L x and T x , while the 
E7 models show numerous small oscillations in L x and T x , that 
increase in amplitude with time increasing (Fig#. The cold gas 
disc is much less extended: its radius is < 1 kpc for k = 0.1, and 
<10 kpc for k = 1 (slightly larger for flatter shape, and for lower 
7 sf). 

3.3 Mass exchange between the cold disc and the new stars 

The second main result of the new simulations is that most of the 
cold gas disc is consumed by star formation, both for pSF and aSF. 
This result, of general validity, is shown for the particular models 
of Fig. [5]in Fig. [8] where the cold gas mass M c of the noSF models 
(black solid line) is at the end very close to that of the new stars M" ew 
(dashed lines). For example, the E07j®° model in the noSF case has 
a cold disc extending out to a 2.5 kpc radius, of M c = 1.6 x 1O 9 M 0 , 
while the corresponding pSF model at the end has no cold disc, 
and a mass in new stars of M“ ew ^ 1.7 x 10 9 Mq (a bit more gas 
has cooled in the pSF case than in the noSF case). Similarly, this 
same galaxy model, with aSF and t/sf = 0.1, at the end has a cold 
disc of just M c = 5.7 x 10 7 M o , and Mf w = 3.8 x 1O 9 M 0 . This 
M" ew is larger (~ double) than for the pSF case, a fact that is partly 
accounted for by the injected mass from the new stars in the aSF 
models and mostly by the delay in the final cooling episode after 
which a galactic wind is established (Fig. [5j; this delay keeps the 
galaxy gas-rich and starforming for a longer time (which is proved 
also by the lower escaped mass M esc in the aSF than in the pSF case, 
see Tab#. 

In general, the LM models with SF show a final M c that can 
range from 0 to 10 9 Mq, but is typically < 10 8 Mq (Fig. [5] and 
Tab. [3]). The final mass in the newly formed stellar disc is of the 
order of M" cw =a fewxlO 9 Mq; thus, the M" cw values reach almost 
those of the M c of the corresponding noSF models (Tab. 0 Fig# 
temperature below 2xl0 4 K. 

Also in the HM models, similarly to what found for the LM 
models, M" ew at the end of the simulations is of the order of the M c 
values of the noSF models. A larger massive disc remains in the 
HM models, though, even when including SF (e.g.. Fig. [8] an d0- 
Their final M c is = 10 9 Mq, for t / s f = 0.1 (= 1/50 of M c without 
SF), and (8 - 12) x 10 9 Mq, if r] SF = 0.01 ( M c = 1/5 of M c without 

2 We recall that, in the aSF scheme, the mass injected back into the ISM 
by the new stars can itself condense and produce new stars, that in turn 
contribute to the computation of M* ew . 


© 2015 RAS, MNRAS POQ.ftp)] 





10 A. Negri, S. Pellegrini, (PL. Ciotti 



5 10 15 20 25 30 


5 10 15 20 25 30 


5 10 15 20 25 30 


R (kpc) 


R (kpc) 


R (kpc) 


<1 — 

I (r b 10“ 4 10“ z 10° 10 2 


Figure 6. Meridional sections of the density distribution of the new stars formed up to 13 Gyr (in units of Mq pc for the same models of Figs. [5]and[8] 
EO7jg 0 (left panels), F07™ (middle panels), and £07^?° models (right panels). Top and bottom panels refer respectively to pSF and aSF simulations, with 
n SF = o.l. 
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Figure 7. Time evolution of the X-ray luminosity Lx and X-ray emission weighted temperature T\ for the EO4j| 0 and E07j-|° models (top panels), and for the 
E04“® ! and E07 ^q j models (bottom panels). For the aSF models, //sr = 0.1. See Sect. 


3.2 


for more details. 
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Figure 8. Left panels: time evolution of the cold gas mass M c (solid lines), and of the cumulative mass in newly formed stars M" ew (dashed lines). Right 
panels: time evolution of the SFR rates. The models are the same of Fig. [5] The black line refers to the noSF model, green and red lines refer to pSF and aSF, 
respectively, with )/sf = 0.1. Note how M c is much reduced in the SF models (green and red solid lines) with respect to the noSF models (black lines), and how 
at the same time M) lcw (dashed lines) reaches close to the noSF M c values. See Sect.|3.3|for more details. 


SF), independently of pSF or aSF. As found for the LM galaxies, in 
the aSF models A/" ew is larger than in the pSF ones, and can be even 
larger than M c of the noSF models (Figs.[8]and[9]l. This time, though, 
the difference is lower, of - 20%, and corresponds to secondary SF 
sustained by the material injected by the new stars. The final mass 


in the newly formed stellar disc is of the order of M" ew =a fewx 10 10 
Mq (Tab.[3|. Flowever, it must be noted that the structural properties 
of the HM models may not be very realistic: ETGs of this mass with 
the isotropic rotator kinematics are quite unlikely, since they are not 
observed to rotate this much (jEmsellem et al.|201 1 j>. Real massive 
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t]sf = 0 ?7sf = 0.01 t]sf = 0.1 



Figure 9. Top panels: final values of the cold gas mass, for the models in Tabs. |2|5| without SF (left panel), and with aSF, for the two different ?7 sf values 
adopted here. Middle and bottom panels: final values of the stellar mass in newly formed stars M" ew , of the SFR normalized to the original M*, and of the mean 
formation time of the new stars (?)* ew , calculated from the beginning of the simulation, for the same models with aSF in the top panels. All values are given in 
Tabs. |33 See Sects. l3.3l and l3Al for more details. 


galaxies will have less rotation than adopted here with k = 1, and 
then they will have less cold gas and less new stars. 

More realistic massive and rotating ETGs are those of Sect. |3.2| 
with aSF. For these, if k = 1, the final cold gas mass is again largely 
reduced when introducing SF, by a factor of =; 50 (E4) and =; 100 
(E7), when 7/ SF = 0.1. At 13 Gyr, is fewxl0 9 Mg for r/ SF = 0.01, 
and fewxl0 s M o for t^sf = 0.1; Mf w is = 2 x 10 10 M o , for both 
77 sf values. If k = 0.1, the final M c is reduced even more than for 
the k = 1 models (Fig. [9]l, for both the E4 and E7 shapes: M c is a 
few x10 8 Mq (jisf = 0.01), and < 10 8 Mq (t/sf = 0.1). Note that the 
final cold mass M c , without SF, is similarly large for the k = I and 
k = 0.1 models (Tabs.[4]an d[5j Fig . [9]), but it is instead lower for the 


k = 0.1 models than for the ( = 1 ones, when adding SF; one could 
conclude that in the ( = 0.1 models SF is more efficient, because it 
takes place in a cold disc which is smaller and denser, or that the 
evolution is different, if less gas overall cooled (the heating was 
more efficient). The second hypothesis is supported by the finding 
that, at the end of the simulations, M” ew is similar for k = 1 and 
k = 0.1 (e.g., 2.4 and 2.3 x 10 10 Mq for the E4 models, respectively 
with ( = 1 and k = 0.1, for t/sf = 0.1), and similarly larger than M c 
of the corresponding noSF models (e.g., 2.0 and 1.9 x 10 io Mq, for 
the E4 models respectively with k = I and k = 0.1). 

In conclusion, SF seems to be an important mechanism to solve 
the worrysome feature of massive cold gas discs in rotating ETGs 
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without SF. In addition, we found the following trends, at the end 
of the simulations (see also the top two rows of panels of Fig. [9}: 
M c increases with M„ and has a large spread for LM models, that 
increases for larger r \ SF ; MJM, is roughly constant with M„, with 
a spread reaching down to values much lower than this constant 
for LM models; at any M„ M c decreases when t/sF increases. M” ew 
and M" cw /M» increase with M„ with a spread for LM models, and 
are roughly independent of rjsF- This latter result means that more 
massive ETGs have been overall more efficient in forming stars 
from recycled material, considering their past = 11 Gyr. 

3.4 The SF history 

SF is not only responsible for the final mass budget of M c and M” ew , 
but also for the rate of SF (SFR), its evolution, and the (fiducial) 
age of M” cw ; that is, for the SF history that we discuss here. Star 
formation shows a different evolution in LM and HM galaxies, in 
pace with their different ISM evolution (see Fig. [8}. In LM galax¬ 
ies, where major cooling episodes take place recursively, the SFR 
peaks suddenly in correspondence of these cooling events due to 
the increased density of cold gas, and declines right after, going 
eventually to practically zero either if the cold gas reservoir is com¬ 
pletely consumed or the cooling episode is followed by a global 
wind (Sect. ED- Thus, in general LM galaxies have a larger and 
peaked SFR in their past (reaching 1-2 Mg yr -1 ), when the rate of 
stellar mass losses due to the original stellar population was much 
larger, and a low SFR (< O.IMq yr -1 ) at the present epoch. Typical 
values at 13 Gyr can be as low as zero, and at most as high as 0.45 
Mg yr -1 (Tab. [5}. These values compare well with the current rates 
estimated for the ATLAS 30 sample, that range from « 0.01 and 3 
Mg yr , with a median value of 0.15 Mg yr -1 for ETGs of a e g 
comparable to those of the LM models run here i jDavis et al.|2014fr . 
HM galaxies, instead, show a more regular and steady production 
of cold gas, and so is their SFR (Fig. [8}; at 13 Gyr, their SFR is 2-3 
Mg yr -1 (Tab.[2] but we recall these models are not realistic). In the 
intermediate mass models (Sect. [X2| l, the SFR is larger for the E4 
models without oscillations in Lx than for the E7 ones. Overall, the 
SFR at 13 Gyr is similar for k = 1 and A: = 0.1: for k = 1, the SFR 
ranges from 0.42 to 1.1 Mg yr -1 (i.e., it is larger than for the LM 
models of same k = 1); when k = 0.1, the SFR ranges from 0.45 to 
1.0 Mg yr -1 . Figure[9]summarizes the present epoch SFR behaviour, 
as a function of M„ for the aSF models. To compare models of 
different mass, the plotted quantity is the SFR normalized to M,. 
This quantity is ~ (1 - 4) x 10 -3 Gyr -1 , quite independent of r/ SF , 
for the intermediate mass and HM models; it can be largely varying 
for the LM models (from 4 x 10 -4 to 6 x 10 -3 Gyr -1 , for t/sf = 0.1. 
and from 8 x 10 -5 to ~ 10 -3 Gyr -1 , for = 0.01). On average, the 
present epoch SFR/AT, increases with M„ a result similar to that of 
the previous Sect. |3.3| where we found that the integration over time 
of the SFR, that is M" ew , and increase with M». 

Tables[2]and[3]also list the mean formation time (f)" ew of the 
new stars since the beginning of the simulation (to = 2 Gyr), and the 
mean star formation rate (defined as M/ ew /(t)" ew ), both at 13 Gyr. 
(r)f w as a function of time is defined as: 

(t)T w (t) =--- f (t' - t 0 ) SFR(t') dt', (9) 

Mf w (f) J , 0 

where the SFR(f) is the volume integrated, instantaneous star for¬ 
mation rate. The values of (f)” ew at the present epoch are shown in 
Fig.[9]for the aSF cases. They range from 2.4 to 6.6 Gyr, for the LM 
models, with a tendency to be lower for the aSF than for the pSF. 
and for the larger t] S p. The HM and intermediate mass models have 


(f)" ew in a narrower range (from 4.2 to 5.5 Gyr), independent of the 
aSF or pSF, and are again lower for the larger 77 SF (Fig.|9j. Therefore, 
the new stars in LM models have a larger spread of formation times, 
and correspondingly a larger range of ages at the current epoch, than 
in larger mass models. 

For the aSF, for rj sf = 0.01, the upper envelope of the more 
recent formation times (largest (t)" ew values) is - 5.5 Gyr, roughly 
independent of M„ while the lower envelope (lowest (f)"™ values) 
increases with M,. A similar trend is shown by the HM and interme¬ 
diate mass models, for = 0.1, but shifted towards lower (t)” ew 
values (< 4.5 Gyr); the spread of (f)” ew for LM models keeps instead 
large, extending to even lower (t)" ew values of = 3 Gyr. The decrease 
°f (t)” ew for larger p SF , i.e., for a more efficient SF, obtained here 
from hydrodynamical simulations provides support to the similar 
conclusion drawn in galaxy formation models focussed on chemical 
evolution: these assume that the efficiency of star formation is an 
increasing function of mass, which for them has - among others 
- the consequence that smaller galaxies continue to form stars for 
longer periods, while more massive galaxies host an older stellar 
population (e.g., |Matteucci|l994) >. 

It is interesting that (f)“ ew decreases slightly, or remains similar, 
when switching from pSF to aSF; in fact, one could expect a more 
extended SF, due to the injection of material from the new stars, for 
the aSF. Evidently, the net results of the injection of both mass and 
energy is to produce the bulk of SF at roughly the same time, for the 
pSF and aSF schemes; only, the final mass in new stars is larger for 
aSF. 

The mean age of the newly formed stars in Gyr is 13-to-(l)” ew - 
This age is in any case larger than ~ 5.5 Gyr, for the aSF. 


3.5 The SF recipe and the Kennicutt-Schmidt relation 

From eq. Hit follows that SF is favoured in regions of high density, 
and relatively short t sF . These conditions are naturally present in low 
temperature regions; accordingly, the simulations show that most of 
the SF takes place where the gas temperature is < 2 x 10 4 K, i.e., in 
the cold rotating disc. 

Quite obviously, our results depend somehow on the adopted 
prescription for the SF, while a realistic, physically consistent, de¬ 
scription of SF is not available. Yet, it is interesting to ask at what 
level our implementation of SF recovers some basic observational 
features of SF. More than half a century ago |Schmidt| ( |1959l ) con¬ 
jectured that the rate of SF should vary as a power law of the gas 
density, and later |Kennicutt|{T998fr suggested the parameterization 
in terms of mass surface densities Zsfr = AL N ., , valid for star- 

forming galaxies, with N = 1.4 - 1.5 ( |Kennicutt & Evans[201 2) . 
We investigate here what scaling (if any) the adopted SF recipe of 
eq.|T]produces in terms of Esfr vs. X gas , when implemented in the 
simulations. 

Figure [T0| shows for a few representative models at the end of 
the simulation the relationship between Zsfr and Z gas , calculated 
separately for the hotter (T> 2 x 10 4 K) and colder (T< 2 x 10 4 K) 
phases of the gas; E S fr is the face-on projection of p SF . The figure 
shows for comparison a recent scaling for the Kennicutt-Schmidt 
relation, derived by Davis et al. (2014) for the ATLAS 30 sample of 
molecular gas-rich ETGs (this has a slope of N = 1.4). Figure [T0| 
tells that SF in the gas takes place basically in two regimes, and 
by far most of it is due to cold gas. The behaviour of the cold gas 
shows two remarkable features: the first is that a scaling similar to 
that of the Kennicutt-Schmidt relation is followed quite closely by 
the blue (cold gas) points; the second is that even the normaliza¬ 
tion is reasonably well reproduced by our scheme of “active SF”, 
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Figure 10. Projected SFR versus mass surface density of the gas (at T < 2 X 10 4 K, for the blue symbols, and at 7 > 2 X 10 4 K, for the red ones), for two aSF 
models: E07 3 * 0 on the left, and EO7~| 0 on the right. Full circles refer to //si = 0.01 and crosses to //sr = 0.1. The lower t ;sf = 0.01 coefficient for the SF recipe 
(eq.fl) reproduces better the empirical Kennicutt-Schmidt scaling, that is shown by the solid line (with a slope of N = 1.4, and normalized as in [Davis ct al.| 
1 20141 for their ATLAS 30 sample of molecular gas-rich ETGs). The blue dotted and dashed lines show the best fits for the cold gas only (coefficients in the 
frames in the lower right comers). See Sect.|3.5|for more details. 


when j/sf = 0.01 (models with tjsf = 0.1 correspond instead to a 
normalization that is too large). In particular, the cold gas in the 
logZ S F R -logE gas plane follows a trend best-fitted by a line of slope 
of a; 1.4 - 1.5, a fact that can be explained as follows. SF in the 
cold gas is most likely regulated by the dynamical timescale t dyn , 
being f coo i very short at such temperatures and densities; thus, one 
roughly expects that /)sf oc p L5 . A slope for the XsFR-^gas relation 
close to the empirical one then follows when neglecting the differ¬ 
ence between volume and surface density, which is not unreasonable 
considering that the cold disc resides mainly in a thin region above 
and below the equatorial plane, extending vertically for just a couple 
of gridpoints (< 200 pc; see also Kennicutt & Evans|2012| l. For the 
hotter gas, the red points in Fig.[l0]are invariably fitted by a steeper 
line, with slope = 2; this is explained by a longer cooling time than 
the dynamical time, and then by psF « p 2 . 

Overall, Fig.[l0]proves that the recipe of eq.[T]is very reason¬ 
able for SF, and, together with the chosen range of normalization 
for rj SF , it provides a scaling close to that observed for normal to 
starburst galaxies. These considerations add strength to the overall 
results obtained in this study. 


4 CONCLUSIONS 

In this work we have explored the effects of SF on the evolution of 
rotating hot gas flows in early-type galaxies; our previous (N14) high 
resolution 2D hydrodynamical simulations of such flows, run for 
axisymmetric two-component models, were successful in producing 
the observed L x and jT x of flat and rotating ETGs, but also revealed 
the formation of very massive cooled gas discs. To study the SF 
effects, we performed hydrodynamical simulations using the same 
numerical code of N14, where SF with a Salpeter IMF is inserted, 
following a simple recipe depending on the gas density, and the 
cooling and dynamical timescales of the ISM (eq. [TJ, with two 
possible values for the efficiency of SF (t 7 SF = 0.01 and 0.1). We 
considered the new stellar generations to be passive (a pure sink of 
gas), or active (the new stars contribute mass and energy to the gas 
flow, via stellar winds and SNIIs). The new simulations have been 
run for a subsample of the N14 models, of high, intermediate and 


low galaxy masses, of E4 and E7 shapes, with the ordered stellar 
rotation described by the isotropic rotator, or with the rotational 
velocities reduced by a factor of ten. 

It is found that subsequent generations of stars are formed from 
the cold gas that accumulates in the equatorial plane, and that most 
of the extended and massive cold disc found by N14 is consumed 
by this process. In particular, the main conclusions are summarized 
as follows. 

• Remarkably, both for the passive and the active SF implemen¬ 
tation, we confirm the results of the previous investigation without 
SF, concerning the trends of L x and T x of the hot ISM with galactic 
rotation. As found by N14, L x is lowered by rotation: in low mass 
ETGs this happens because rotation favours the establishment of 
global winds, and in medium-to-high mass ETGs because rotation 
lowers the hot gas density in the central galactic region. The average 
temperature T x is also lower in intermediate/high mass ETGs, when 
rotation is important, due to the lower contribution of the central 
regions (usually hotter and denser in non-rotating systems), while it 
can be higher in low mass ETGs, if rotation triggers a wind, due to 
the thermalization of meridional winds. 

• The robustness of the N14 results even after the addition 
of SF is explained differently for HM and LM galaxies. In HM 
galaxies the evolution of the hot ISM is only marginally sensitive 
to the presence of SF, and L x and T x evolve smoothly as without 
SF: a massive cold disc forms early, it is mostly consumed by SF. 
and around it the hot ISM evolution proceeds almost unaltered by 
SF (even in presence of energy and mass injection by SNII’s). In 
less massive models, instead, the evolution of L x is not smooth, but 
shows many peaks corresponding to major cooling episodes, as in 
N14. The time-occurrence of the peaks depends on SF, thus SF can 
induce variations of L x at a certain chosen epoch; however, these 
keep within the already large range of L x values typical of these 
masses, even without SF, being the gas flow very sensitive to many 
factors (as already found previously, e.g., |Ciotti etal.|199lj N14). 

• For what concerns the global gas mass budget, at the end 
of the simulations the cold gas mass left in the equatorial disc M c 
on average increases with M„ with a larger spread for LM models; 
at fixed M„ M c is lower for lower stellar ordered rotation, and for 
larger In addition, the ratio MJM, is roughly constant with M„ 
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with a spread reaching down to values much lower than this constant 
for LM models. In any case, M c < 2 x 10 9 Mq for galaxies with 
M. 4 x 10 u Mg. Typical values for M c are <a few xlO 7 Mq for 
LM galaxies (with the maximum rotational level of k = 1), with 
exceptions of M c - 10 9 Mq; and a few xlO 7 Mq (if k = 0.1) to 
3 x 10 9 Mq (if k = 1) for ETGs of intermediate mass. These values 
compare well with those recently observed I Young et al. 2014| |Seira| 
|et al.|2 012 20141. In particular, they can explain the observation 
that massive, fast-rotating ETGs often have kinematically aligned 
gas, independent of environment ( jDavis et al.|2011 [|2013| >. 

• The mass in newly formed stars M" cw , and the ratio M" ew /M,, 
increase with M„, again with a large spread for LM models, and are 
roughly independent of rj sf- This result means that more massive 
(rotating) ETGs have been overall more efficient in forming stars 
via recycling of their stellar mass losses, during the past ^ 10 Gyr. 
The mass in secondary generations of stars is (1 - 6) x 10 9 Mq (if 
k = 1) for LM models, and ^ 2 x 10 10 Mq (for both k = 1 and 
k = 0.1) for intermediate mass models. These should reside mostly 
in a disc, as most fast rotator ETGs possess, and be related to the 
birth of a younger, more metal rich disky stellar component that 
is indeed observed l |Krajnovic et al.|2008[ (Cappellari et al. 2013[ >. 
They should not be recognized as a “young” stellar population in an 
absolute sense, though, since most of M" cw formed a few Gyr ago, 
as can be evaluated from the evolution of their SFR. 

• The time evolution of the SFR depends on the mass of the 
galaxy. LM galaxies (M, = 1.25 x 10 11 Mq) have a larger and 
peaked SFR in their far past (reaching 1-2 Mq yr~'), when the 
rate of stellar mass losses due to the original stellar population was 
much larger, and when major cooling episodes were frequent. In 
fact, the average age of the new stars of LM models ranges from 5 
to 8 Gyr (for the aSF scheme). At the present epoch, their SFR is 
low (< 0.1 Mq yr~* typically), with the full range of values going 
from zero to ^ 0.7 Mq yr _1 . These results agree nicely with the 
low degree of SF and young stellar populations that is detected only 
in fast rotators, in the ATLAS 30 sample l Kuntschner et al.|2010 


|Sarzi et al.|20 13 1; also, they compare well with t he current estimated 
rates, whose median value is » 0.15 Mq yr 1 i Davis et al,|2014 1. 
More massive galaxies, instead, show a more regular and steady 
production of cold gas, and so is their SFR; the average age of their 
new stars ranges from 5.5 to 7.0 Gyr (in the aSF scheme). At 13 Gyr 
the SFR is ^ (0.4 - 1) Mq yr -1 , larger than for LM models, and 
somewhat larger than observed (e.g., [McDermid et al.|2015| ; this 
could be fixed by assuming a lower or a lower ordered stellar 
rotation, at larger galaxy masses (Fig.[9]i. Thus, our models may still 
be rotating too much and producing too much cold gas, at galaxy 
masses > 2 x 10 11 Mq. The current SFR/M*, instead, remains quite 
constant from LM to intermediate mass models, with a larger spread 
for LM models (as already found for other properties). 

Finally, the SF recipe adopted for this work proved to be a 
reasonable one, given that it reproduces the slope of the Kennicutt- 
Schmidt relation, and even the normalization if j/ S f - 0.01. 

Overall, we can suggest an origin (mostly) in the SF from 
cooling hot gas, for the presence of cold gas phases kinematically 
aligned with the stars, and for the low-level degree of SF, all features 
detected only in fast rotators in the ATLAS 30 sample. 

In the perspective of the galaxy evolution, we finally recall that 
the precise knowledge of the amount of gas flowing towards the 
galactic centre is of great importance for a proper study of feedback 
effects from the central black hole in rotating galaxies (e.g., [Novak] 
|et al.|20lT[ [Gan et al. 20141. A major implication of the present work 
is the fact that, if the gas produced by stellar evolution forms new 
stars, the amount of gas in principle available for fuelling the central 


supermassive black hole is much lowered. Of course, such fuelling 
would be impossible in presence of rotation and absence of viscosity. 
Disc instabilities, though, could break axial symmetry, and allow 
for the gas infall onto the black hole, possibly escaping SF. Phe¬ 
nomenologically, the effect of gravitational instabilities in the disc 
can be modelled as an effective gravitational viscosity (e.g., |Bertin &| 
|Lodato|2001||Hopkins & Quataert|20 1 I [|Rafikovj20 1 5| , that favours 
accretion of cold gas towards the centre. Indeed, viscosity (that 
could be due for example to MRI) is another ingredient not present 
in the current simulations, but that could affect the disc evolution. As 
the code is axisymmetric, we cannot follow the complex physics of 
(non-axisymmetric) disc instabilities that could be present. Although 
self-gravity of the gas is not considered in our simulations, yet for 
some models at a selection of times we computed the fiducial value 
of the radial profile of the stability parameter Q(R) = c s k R lnGT. z , 
where k R is the local epicyclic frequency of the galaxy, and c s and 
2 C are respectively the vertically mass-averaged values of the ISM 
sound velocity, and the surface density of the cold rotating gas. The 
Q values were computed for a layer Az - 200 pc thick above and 
below the galaxy equatorial plane; these Q turned out to be quite 
independent of the adopted thickness value Az, being dominated by 
the cold disc. These fiducial Q values maintained invariably larger 
than unity (indicating stability) when SF is allowed, while Q < 1 on 
a central ^kpc-size region, in absence of SF. From this preliminary 
and qualitative analysis, it follows that SF and (^-instability are 
both present in the gaseous rotating discs. The relevant question is 
then if and what of the two processes is dominant in the disc. In a 
future work we plan to extend further the present investigation by 
implementing in the code the cooperating effects of SF and mass 
discharge from the disc to the centre. In this new study, we will also 
consider the neglected late-time mass return from the newly born 
stellar population, that is the mass injected by stars of mass < 8 Mq; 
this input could be modelled with the method described in [Caluraj 
[etaL| ( (20T4) . 
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Table 2: Relevant quantities at t = 13Gyr for the rotating models (k = 1) with cr e8 =300 km s' 
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Notes. (1) Model name. (2) Star formation scheme adopted (x=noSF, P=pSF, A=aSF). (3) Star formation efficiency. (4) - (5) Cumulative mass injected by the evolution of the original stellar 
population, and mass escaped from the numerical grid, respectively. Differences in Mj n j for models of same Lg are accounted for by different sampling of the stellar density profile over the 
numerical grid. (6) Total ISM mass retained within the galaxy (7) - (8) ISM mass with T > 10 6 K, and T < 2 x 10 4 K, respectively. (9) Mass of the new stars. (10) Cumulative mass ejected by 
type II SNe. (11) — (13) Mean formation time of the new stars, mean SFR, and SFR at the end of the simulations. (14) ISM X-ray luminosity in the 0.3-8 keV band. (15) - (16): SNIa and SNII 
kinetic energy input per unit time (Sect. 2.2.2). (17) ISM X-ray emission weighted temperature in the 0.3-8 keV band. 











Table 3: Relevant quantities at t = 13Gyr for the rotating models (k = 1) with cr e8 =200 km s' 
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